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Abstract 



^ \ We discuss direct CP violation in three-body decays of B meson 

such as B — > K ± 7t + tt~ , which involves various intermediate resonance 
states, such as B^ — ► K*^ and B ± — > pK ± , where K* represents 
the ttK resonance. Due to the large final state interaction phases of 
the resonances, the CP asymmetry can be as large as the 25% level 
near the kinematical region where the K* and p resonances overlap. By 
examining a Dalitz plot of this mode and combining a measurement of 
the branching ratio for B^ — > p K, it is possible to extract the weak 
phase, (p3=axg(—V u dV* b /(y c( iV* b )), of the Kobayashi-Maskawa matrix. 
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The origin of CP violation is one of the unresolved problems, although CP 
violation has been observed in the K°-K mixing system for 30 years. CP vio- 
lation occurs due to the complex phase of the Kobayashi-Maskawa (KM) matrix 
in the standard model (SM). The main purpose of B- factory experiments is 
to measure the phases of the KM matrix, defined by (fii=axg(—V cd V* b / (V td V t l)) , 
(j) 2 =arg(-V ud V* b /(V td V t l)) and (j) 3 =arg(-V ud V* b /(V cd V; b )), and to check the unitarity 
nature of the KM matrix || Q . It is known that <\>\ can be measured in a so-called 
"gold plated" B° — > J/ipK s mode Q. This mode is expected to determine 4>i to 
good accuracy. On the other hand, measurements of 02 and 03 are considered to be 
rather difficult from the experimental point of view. Various methods to measure the 
weak phases are discussed in the literature 0. In order to observe the CP- violating 
effects it is necessary to have different CP-conserving and CP-violating phases in 
an amplitude. The CP-violating phase is provided by the phase of the KM matrix. 
For the neutral B meson the CP-conserving phase is provided by the phase of B°B 
oscillation, whereas for the charged B meson the CP-conserving phase is dominantly 
supplied by the final- state-interaction (FSI) phase. In general, it is difficult to reli- 
ably evaluate the FSI phase because of hadron dynamics. However, we can reliably 
calculate the strong phase in the case of resonances by using the Breit-Wigner form. 
In addition, the resonance provides a large CP-conserving phase, which may lead to 
a large direct CP- violation effect || and the resonance enhancement of CP violation 
is discussed in several decay modes J7|, ||. 

In this paper we propose another method to measure the angle 03 using the 
direct CP violation in three-body decays of the B meson, such as the i? ± — > K ± 7t~tt + 
decay mode. Several resonances contribute to this decay mode, the effects which 
might enhance direct CP violation. There are two channels of two-body intermediate 
states in this decay mode. One is B^ — > K^p — ► K ± 7t~7t + ; the other is B^ —>■ 
K*^ — > K ir~ir + , where K* denotes the ttK resonances summarized in Table.|l] ||. 
The B^ — > K^p process involves equally tree (T) and QCD penguin (P) diagrams, 
whereas B ± — > K*^ has a dominant contribution from the QCD penguin (P') 
diagram. § The amplitudes can be written by 
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A(B~ - <0 = Pi i 



§In principle there may be a tree-diagram contribution to B — * K*w through a rescattering 
process, B ± — ► K^p — ► K*-k^ and this contribution will be discussed later. 



b p = A(B~ -»• K'p) = Te iST e- i(t>3 + P e i5p , 



where St, 8p and S' P . represent the strong phases, and the tree-diagram contribu- 
tion has the weak phase 03. Concerning the corresponding amplitude for charge- 
conjugated modes al and b p , only the sign of 03 should be changed. These amplitudes 
are normalized as 
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BR(B~ -> K~p) = \b p 
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Although these modes have not been observed experimentally, the CLEO collabo- 



ration reported the branching ratios for the B —* ixK modes as follows [jTU 

= (2.3±^±0.2±0.2) x 1(T 5 , 
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This result indicates that the T, P, P' amplitudes for the B + — > K tt + and B° — > 
K + 7r~ processes are the same orders of magnitudes |TT|| : Since the B° — > tt + K~ and 
B~ — ► tt~K transitions are essentially the same as the B~ — > K~ p and B~ — > K*7r~ 
transitions, respectively, these intermediate contributions can significantly interfere 
with each other in the B^ — > K ± n + ir~ decay modes. The intermediate resonances, 
which subsequently decay to two pseudo-scalars, have large FSI phases. We assume 
that these resonance contributions can be treated by the Breit-Wigner form: 
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s' = (p n + +p n -) 2 and the m^s and IYs are the masses and 



s — m% 

where s = (p K - + £>7r+) 2 , 

widths of the resonances, respectively. 

Assuming that the B~ — ^ K~7t~tt + decay amplitude is dominated by the above 
resonances, we can write the decay amplitude in a narrow-width approximation as 
follows: 
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where B^ 71 " and B^* denote the branching ratios for p — > 7T7t and if* — ► Kir, respec- 
tively, and (31 is defined by 



1 - (mj + m k ) 2 /m 2 )(l - (rrij - m k ) 2 /m 2 ) . (8) 

The function G describes the angular component of the amplitude for the resonance 



decay || [12|]. The explicit form depends on the spin of the resonance, and is given 
by 

| \[\, spin 0, 

&(z) = I s[\z, spin 1, (9) 

^ yj\(Zz 2 - 1), spin 2, 

where z(z') denotes cos#(cos#'), and 6(9') is the angle between the momentum of 
the B and K~(ir~) mesons in the center-of-mass frame of the 7r + 7r~ (tt + K~) pairs. 
z(z') can be written as a function of s, s' as follows: 

2ss' + s 2 — s(m 2 B + m 2 K + 2ml) + (m 2 B — ml)(m 2 K — ml) 



^J(s - (m\ + m n ) 2 )(s - (m K - m 7r ) 2 )(s - (m B + m 7T ) 2 )(s - (m B - m n ) 2 ) 
\/s 7 (2s + s' - m B - m\ - 2m 2 ) 



(10) 



(s' — Am 2 )(s' — (tub + r n^K) 2 )(s' — (tub — thk) 2 j 
For simplicity, we hereafter consider the case that only the if 2 (1430) and p res- 
onances contribute to each irK and nir channel, respectively' We later discuss more 
general formula which includes all of the intermediate resonances. The differential 
decay width for B~ — > K~tt~tt + can be written as 

dT(B- -»■ K-TT-n + )/dsds' (12) 

2 
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Re(f( K ^\ s ,s')) + \b p \ 2 Re(g^(s,s')] 

+ Re (a*. 6;) Re (h^' p \s, s')) - Im (a K * b* p ) Im (h^' p \s, s')) . 

The distribution on a Dalitz plot is described by the functions /, g and h in this 
formula, and are written as 

2t)KttT}Ktt 
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^The if|(1430) resonance overlaps with the Kq (1430) resonance but we neglect it here. We will 
discuss this contribution later. 
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The above expressions are defined for general meson states, fi—p, /o, or / 2 , which 
decay into the 27r state since we will later discuss the B^ — > /i7r ± — > 7r + 7r~i^ mode 
in addition to the 5 ± — > pK^ — ► 7r + 7T~K ± mode. We can obtain the differential 
width for the charge-conjugated mode by simply changing the amplitude a,i(b p ) into 
Hiipp). In Fig.|I](a)-(d) we show the kinematical functions Re(f( K %' K ^), Re(g ( - P ' p ^), 
Re(/?,(- K 3'P)) and Im(/i (if 2>p)) i n the region R = {0 < s, s' < 3GeV 2 }. It can be seen 
that the signs ofKe(f^' K ^') and Re(g < - P ' p )) are always positive, whereas Re(/i^2>P)) 
and Im(/i^2'W) change their signs around the K^l^O) and p poles. Note that the 
i^2*(1430) and p contributions can overlap around the region s ~ m 2 K *, s' ~ m 2 and 
that there is a large interference in this region. 

Let us discuss how we can determine the parameters in Eq.(|l2|). Since one of 
the strong phases is irrelevant, we have six unknown parameters (T, P, P', 5t — S' P , 
5p — S'p and fa) in this formula. We can determine these parameters as follows. In 
the kinematical region s ~ m 2 K * and s' 3> m 2 the K\ contribution dominates in the 
decay width. We can extract the amplitudes \q,k*\ and \Uk*\ by retaining only the 
first term in the Eq.(|12|) and the parameter P' can be determined as P' = \a,K*\- 
By examining the distribution on the Dalitz plot in the overlapping region s ~ m 
and s' 
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m~ 



we can extract Re(a^-*6*), Im(a^*6*), Re(a^*6 ) and Iia( a K*b). Since 

At A " At A " 

we have already known |oj<-*| and |a^*|, we can determine Ke(b p e lSp ), Im(b p e tSp ), 
Ke(b p e lS p) and Im(b p e t5 p). There are five unknown parameters in b and b but we 
have only four constraints. If we combine the branching ratio for B~ — > p~K, whose 
amplitude is written as 



A(B~ -»■ p~K) 



■V2Pe iSp , 



(16) 



by using the isospin symmetry |T3| , we can therefore fix all the parameters, including 
the weak phase, 03. Note that we assume that there is no tree-diagram contribution 
to the B~ — > K i\~ amplitude in Eq. ([!]). In general the tree-diagram contribution 
can be induced through the rescattering process such as B — » pK~ — > K 2 7f~ ■ If we 



take into account the tree-diagram contribution, which can be written as T' e %&T e~~ l ^ z , 
we have two unknown parameters, 5' T and T", in addition to the above six unknown 
parameters. We cannot determine 03 since there are only seven independent con- 
straints even if we combine the branching ratio for B ± — > p^K. Although it is 
hard to calculate the corrections from the rescattering effects quantitatively, it is 
estimated to be of order 10 % for the B — > nK modes [14]. Unless this rescattering 
effect dominates in the amplitudes, we can expect to determine 03 with reasonable 
accuracy. 

The CP-violating quantity Acp is defined as 

_ dT(B- ->■ K-7i + n-)/dsds' - dT(B + ->■ K+tt+tt-))/ dads' 
CP = dT(B- -> K-ir+ic-) /dads' + dT(B+ -> K+ir+ic-)) / dads' ' ^ ^ 

where we can write the denominator and numerator as follows: 
dT(B- -> K-TT + TT-)/dsds' + dT(B + ^ K + n + TT-)/dsds' 

2P' 2 Re (/(**.*S>) + 2 (T 2 + P 2 + 2TP cos 5 P cos 3 ) Re (g M ) 



+ 2P' (T cos 5' P cos 03 + P cos(5' P - 6 P )) Re (h^^ 

- 2P'(Tsm5' P cos(j) 3 + Psm(5' P -5p))Im( K h {K 2' p) ) } (18) 

dr(fi _ -► K-n + TT-)/dsds' - dT(B + -► K + n + ir-)/dsds' 

= -2T [2Psin5 P Re (V P ' p) ) + P'sin5 P Re (V*^) 

+ P'cosc^im (/i^'d) ] sin 3 . (19) 

The CP asymmetry appears as an asymmetry of the distribution on the Dalitz plot 
for B + and B~ decays. Let us consider here the integrated asymmetry in the region 
R. When we integrate the functions Ke(h^ K ^' p ') and hn(h^ K ^' p ') in this region, they 
almost vanish because the functions change their signs around the resonance poles, 



J R (dT(B- -> K-n + n-)/dsds' + dT(B + -> K + n + n-)/dsds') 

f [2P' 2 Re (/(^a*-**)) + 2 (T 2 + P 2 + 2TP cos 5 P cos 3 ) Re (g M )] .(20) 






In order to see how the CP asymmetry can be large, it is necessary to evaluate the 
hadronic matrix elements, such as T, P and P', which we cannot calculate reliably. 
Instead of evaluating the hadronic matrix elements based on a specific hadronic 
model |l||, let us make a crude estimation assuming the relation T = P = P' . We 
can obtain the maximal CP asymmetry by choosing the integration of Eq.fll9|). If 
we define 

J R dsds'\Re[h^^{s,s f 
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J R dsds' [Re (/(*a.*a)( 5 , s')) + 2Re (#0>>")( s , s')j\ ' 

^4r and ^4/ correspond to the optimal CP asymmetry in the case of 5p = 0, 
tip = 03 = 7T /2 and <5p = 5' P = 0, 3 = n/2, respectively" We also calculate 
CP asymmetry using other tnr resonances, such as /o and f% summarized in Table.|2]. 
Although the isospin relation cannot be used for /o, fi and therefore we cannot 
determine 03 from these modes, a large direct CP asymmetry may be observable. 
In Table. || we present the values of Ar and Ai not only in case of K% and p, but also 
in case of other resonances. Although the actual CP asymmetry indeed depends on 
the poorly known value such as 03, 8p and S' P , it is reasonable to expect that the 
CP asymmetry can be as large as 25% level at somewhere around the overlapping 
region. 

It is straitforward to write a general formula which includes all of the interme- 
diate resonances as follows: 

dT(B- -► K-TT-TT + )/dsds' (23) 

= £|a«| a Re (/<*?•*?>(*, a*)) +I>| 2 Re(<^( S , S ')) 

i i 

+ 2£{Re(a 4 a *) Re (/^'^ (,,,')) -Im(a 4 a*) Im (/<*W(«, s'))} 

i<j 

+ 2 J2 {Re (h b*)Re(g^\s,s'))-Im(b t b*)lm(g^\s,s'))} 

i<j 

"if Sp is large enough, it is possible to have a large direct CP violation in B ± — ► K^p. Hence 
we may observe the direct CP violation even in the region s 3> vnr K *^ s' ~ m 2 . 
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Since the mass of Kq (1430) is as same as that of .Kj; (1430), these resonances coincide 
at the kinematical region -\/s ~ 1430 MeV. Even if the masses are degenerate, we 
can determine |ok-*|, |ok-*|, Re(a,K*ciK*) and Im(a,K*Q>K*) by analyzing the Dalitz 
plot in the region s ~ m K i and s' ^> m? since the Kq and i^| contributions depend 
differently on s'. Once ax* and ax* are known up to an overall phase, we can extract 
Re(6j), Im(6j), Re(bj) and Im(bj) similarly as in the case of K% alone by examining 
the Dalitz plot in the overlapping region. In case of the p resonance we can again 
obtain another constraint from B~ — > p~K with the isospin symmetry. We can 
therefore determine all the parameters including 03. In case of /o and J2 resonances 
we cannot determine 03 although a large direct CP violation may be observable. 

Here, we discuss the experimental feasibility to detect this CP violation. The 
detector and accelerator are assumed to be similar to the Belle detector in the KEK 
B-factory 0. The decays of B — > fC^(1430)7r an d pK are used. Both branching 
fractions, i.e., BR(B -► K*n)-BR(K* -► K + n-) and BR(B -^ pK)-BR(p -» tt+tt-), 
are set to 10 -5 . The kinematical region of s < 3GeV and s' < 3GeV is about 
12% of the total phase space, considering the experimental acceptance and decay 
angular distributions. The total experimental acceptance is about 3%, including 
the geometrical factor and other cuts, such as to suppress continuum background 
(e + e~ — > (7) —>■ qq). The signal-to-noise ratio is estimated to be ~2. In the 
case of having a 25% direct CP asymmetry, ~ 300 ftT 1 of the data on T(4s) is 
necessary to obtain a 3a asymmetry. A significant improvement can be expected, if 
the kinematical region is reduced to only around the resonance overlapping region. 

In this paper we consider the direct CP violation in three-body decays, such 
as .B — > K ± 7i + ii^ . The B meson can decay to the three-body final state through 
two different channels, K*^ and pK^, and these amplitudes may interfere with 
the same orders of magnitude near the kinematical region where the K*7r and 
pK ± contributions overlap. The resonances K* and p provide a large final-state- 
interaction phase, which is reliably calculated using the Breit-Wigner form. This 
yields a large CP asymmetry around the overlapping region, and we can expect the 
CP asymmetry to be as large as the 25% level in this region. By measuring the 
distribution on the Dalitz plot and combining the measurement of branching ratio 



for B^ — > p^K, it is possible to extract the weak phase, 03. 
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K* 


spin 


mass (MeV) 


width (MeV) 


B(K* -»■ K-n+) 


IT* (892) 


1 


896 


51 


0.67 


K*(U30) 





1429 


287 


0.62 


isT* (1430) 


2 


1432 


109 


0.33 



Table 1: Summary of Kir resonances. 



fi 


spin 


mass (MeV) 


width (MeV) 


B(/ 4 -► 7T + 7T-) 


p(770) 


1 


769 


151 


1 


/o(980) 





980 


40-100 


0.52 


/ 2 (1270) 


2 


1275 


185 


0.57 



Table 2: Summary of im resonances. 
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fi 


A R 


Ai 


K* 


P 


0.13 


0.087 


K* 


/o 


0.15 


0.10 


K* 


h 


0.26 


0.25 


K 


P 


0.20 


0.22 


K 


/o 


0.26 


0.28 


K 


h 


0.26 


0.27 


K 2 


p 


0.26 


0.25 


K 2 


h 


0.29 


0.26 


K 2 


h 


0.26 


0.25 



Table 3: CP asymmetries defined in Eq.([2TD-([22i). Each row corresponds to the 
asymmetries, taking into account only one resonance for each im and ttK channel. 
Here the width for f is assumed to be 70 MeV. 
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pvsK 2 



^0.14 

^0.12 

* 0.1 

0.08 

0.06 

0.04 

0.02 






FIG. 1: Kinematical functions denned in Eq.( 



We take into account the p and 



A^ (1430) resonances for the tttt and Kn channels, respectively, s = (px- + Pw+] 
and s' = (p n - + p n +) 2 . 
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